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Abstract 

We derive Bogomolny equations for an Einstein- Yang-Mills-dilaton-u 
model (EYMD-ff) on a static spacetime, showing that the Einstein equa- 
tions are satisfied if and only if the associated (conformally scaled) three- 
metric is flat. These are precisely the static metrics for which which 
super-covariantly constant spinors exist. We study some general proper- 
ties of these equations and then consider the problem of obtaining axially 
symmetric solutions for the gauge group SU(2). 
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1 Introduction 



The Bogomolny equations have played a central role in the search for solu- 
tions of the Yang-Mills-Higgs (YM H) system of equations in the limit of vanishing 
self-interactions of the Higgs field. Because of their lower order, the Bogomolny 
differential equations prove a major simplification to the problem of finding an- 
alytical solutions for the YMH equations. It is natural therefore to seek similar 
equations in the study of Einstein- Yang-Mills- Higgs and related systems in which 
the coupling to gravity increases the (already high) non-linearities of the cou- 
pled, second order equations. This paper presents several new results along such 
lines. 

The investigation of curved-space Bogomolny equations has interested many 
people. A modified version of the Euclidean Bogomolny equations was considered 
by Comtet |2| very early on while studying Yang-Mills-Higgs (YMH) systems (in 
the Prasad-Sommerfield limit) on fixed, static, curved-space backgrounds utilis- 
ing a spherically symmetric ansatz. Here a system of first order Bogomolny type 
equations was proposed that implied the YMH equations of motion provided 
the metric satisfied a differential constraint. After this, in an important work, 
Comtet, Forgacs and Horvathy || treated the general fixed, static, curved-space 
background with no (spatial) symmetry assumptions. Again first order Bogo- 
molny type equations were found that yielded the YMH equations of motion 
provided the metric satisfied the differential constraint: 

Alnv^ = 0. (1) 

Going beyond Comtet's earlier work, these authors now asked about the com- 
patibility of this constraint with the Einstein equations (here only considering 
the spherically symmetric case). Compatibility was found possible only for a 
very special physical solution: the Higgs field was required to take its vacuum 
value Tr0 2 = 1 and the metric described an extreme Reissner- Nordstrom black 
hole. As we shall see, this first constraint (that Tr0 2 = 1, which may be viewed 
as a sigma model constraint) is fundamental, and appears either implicitly or 
explicitly in most authors work. 

Balakrishna and Wali ||J] , working with an Einstein- Yang-Mills-Higgs (EYMH) 
system with non-minimal RTr<fr 2 coupling, and assuming a conformstatic met- 
ric with flat three-geometry, also obtained the same Bogomolny equations of 
Comtet et al. Here the Einstein equations were found to be satisfied as a conse- 
quence of the Bogomolny equations provided Tr0 2 = 1 and the choice of units 
47rG = 1 was made. The sigma model constraint in fact means that there is 
no difference between minimal and non-minimal gravitational coupling for the 
solutions being considered; the choice of units corresponds to a balance between 
the attractive gravitational and repulsive gauge force of like particles, and will 
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be commented upon further in the sequel. Balakrishna and Wali also examined 
the large and small r asymptotics of the resulting metric, assuming spherical 
symmetry. In a similar vein Cho et ai.[0] considered a (now minimally cou- 
pled) Einstein- Yang-Mills-Higgs-Dilaton (EYMHD) system for a static metric 
with flat three-geometry, their work assuming the constant Tr0 2 sigma model 
constraint. Analogous Bogomolny equations to those found earlier (but incor- 
porating the dilaton) were obtained. Various multipole solutions generalising 
the Majumdar-Papapetrou || [7| electrovac solutions were found, as well as the 
Gross-Perry-Sorkin multi-monopole solution obtained from a five dimensional 
Kaluza-Klein theory [§, g, [TUJ. 



Most recently Forgacs et al. |TTJ have made explicit the recurring sigma model 
constraint noted above by considering the coupling of an Einstein- Yang-Mills 
system to a gauged sigma model (EYM-cr) or an Einstein- Yang-Mills-dilaton 
system with analogous coupling (EYMD-cr). They view this as the infinite mass 
limit of an EYMH system: the infinite mass forcing the length of the Higgs to 
assume its vacuum minimum. The key point is that when coupling to a sigma 
model field n, the differential constraint ([!]) is replaced by 



AlnVkol =Tr(An£> l n), (2) 

and the Yang-Mills and sigma model equations of motion follow from a Bogo- 
molny equation if and only if the metric satisfies this equation. Now (again with 
units 4tcG = 1) they show that ([|) is equivalent to the (00) component of the 
Einstein equations for a general static metric. Although these authors do not 
prove the compatibility of their Bogomolny equation with the remainder of the 
Einstein equations they do present a particular case for which this compatibility 
holds. This particular case corresponds to spherical symmetry and the t'Hooft 
ansatz: the resulting ODE's corresponds with those found in j| noted above, 
which shows that the issue of the compatibility is nontrivial. Most recently Viet 
and Wali fl2| adopting a spherically symmetric ansatz show the compatibility of 
the Bogomolny equation with the equations of motion and Einstein equations of 
a non-minimally coupled EYMH system. Here again a sigma model constraint is 
imposed, and we again encounter the same equations of Forgacs et al. but now 
the full compatibility with the Einstein equations is shown within their ansatz. 

In this paper we follow Forgacs et al. by initially considering an EYM-a 
model on an arbitrary, static spacetime. Our first result is to show that, with the 
assumption of static fields and units so that 47rG = 1, the Bogomolny equations 
(see (0) or (p0|) ) yield all of the field equations if and only if the (appropriately 
conformally scaled) three-metric is flat. This at once answers the question as 
to the compatibility of the Bogomolny equations with the remaining Einstein 
equations and justifies the ansatz of the various authors mentioned above. In- 
terestingly the space-times we find Bogomolny equations for are precisely those 
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static space-times for which super-covariantly constant spinors exist [jnj . This is 
consistent with the arguments of |14] who associate supersymmetric extensions 
to theories exhibiting Bogomolny bounds and suggests a supersymmetric exten- 
sion exists to this theory Having obtained our Bogomolny equations we proceed 
to investigate these in section three. First we show that an auxiliary magneto- 
static problem may be associated with the equations for any gauge group. This is 
particularly relevant for situations where the sigma model field is covariantly con- 
stant and leads to the Majumdar-Papapetrou electrovac solutions noted above. 
To examine the case of non-covariantly constant solutions we focus on the case 
of axially symmetric charge one solutions for the gauge group SU(2); the ansatz 
we adopt readily incorporates the spherically symmetric ansatz of the earlier 
works and we recover these. Here we also give various improved asymptotic so- 
lutions and are able to show that an asymptotically flat, spherically symmetric 
black hole solution with regular gauge fields satisfying the Bogomolny equations 
is necessarily an extreme Reissner-Nordstrom black hole. Section four extends 
our earlier analysis to include a dilaton coupling: up to some physical redefini- 
tions we find that the Bogomolny equations (and consequently their solutions) 
are unaltered. Our penultimate section considers axially symmetric SU(2) so- 
lutions to the Bogomolny equations with charge (or winding number) greater 
than one. Here we can illustrate covariantly constant solutions with arbitrary 
charge but are unable to find solutions to the ansatz of various authors for the 
non-covariantly constant case. We conclude with a brief discussion of our results. 
An appendix comments on the singularities of the solutions under consideration. 



2 The Bogomolny Equations 

Our study of the EYM-a system is based upon the action 



S = Id xy\g 



^—R - - Tr F UU F^ - - Tr D.nD^n + - (Tin 2 - l) 
16nG 4 2 2 V 



• (3) 



Here G is the gravitational constant, R the Ricci scalar associated to the space- 
time metric g^ u and the scalar field n is in the adjoint representation of the 
gauge group with associated field strength F. Indices /i, z/, ... run from to 
3 and we are working with a signature ( — h ++). The Lagrange multiplier A 
of the final term in the action imposes the a-model constraint. This action has 
been considered previously in Ref . []TTJ as the infinite mass limit of spontaneously 
broken gauge theories with adjoint Higgs fields. 
The field equations derived from (|^) are 

1 - T Dj y /\j\F' a ')=[n,irn], (4) 



\9\ 
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D, 



\g\D»n 



Tr D^nD^n j n, 



(5) 



2^ T 



(6) 



where R^ v is the Ricci tensor, is the total energy-momentum tensor associ- 
ated to the gauge field and the cr-model, and T = T 7 ^. (See for details.) 

We restrict our attention to the static, purely magnetic case (Aq = 0) of this 
theory and for the purpose of this paper assume a static metric parameterised 

as 



ds 2 = goo dt 2 + gijdx t dx j = —V 2 dt 2 + ^dx l dx j 



(7) 



Comment on the extension to a stationary metric will be made in the sequel and 
here i,j run from 1 to 3. For such a metric Roi = 0, and we note the identity 
Ej[p 339] 



R 



00 



-VD g V 



-V 3 D h V + V 2 h ij d i Vd j V, 



where O g (h) is the scalar Laplacian with respect to the three-metric gij (or hij). 
Since the time derivative of every field is assumed to be zero we may construct 
a reduced action S, with S — J dtS, and whereQ 



S 



d 3 xVh 



16nG 



(R {3) (h) + 2V* Vi In V - 2 V* In V V, In V) 



XV 



-2 



-V 2 Tr Fi~F ij - - Tr D^n + (Tr n 2 - I) 



(9) 



In this reduced action, and indeed from here onwards, all covariant derivatives 
are constructed from and indices lowered etc. with respect to the three metric 
h^, and R^ 3 ' (h) is the scalar curvature of this metric. We obtain @ by expressing 
R(g) in terms of R^ (h) and the function V. For the reduced action the equations 
of motion for V correspond to the (OO)-Einstein equations (|). 

Our strategy to derive the Bogomolny equations is in the spirit of [0] where 
actions were expressed as sums of squares. First we show that (with an appro- 
priate choice of units) the reduced action may be expressed in the form 



S 



J d 3 xVh ^R (3) (h) -^Trti j v?vf 



+ B 



± 



(10) 



1 With our assumptions of static, purely magnetic, fields we have an energy density Tq 
Cm, these being the final three terms of (ph. 
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where v f is an appropriate combination of the fields and is a surface term 
that will be specified shortly. The key point is that the action is quadratic in v ± 
(and this quadratic is of definite signature say for a semisimple gauge group). 
Thus for any field variation this term of the action has an extremal when = 
for every i and gauge index a. These will be our Bogomolny equations. Further 
analysis will then show the three-metric hij flat. 

Towards establishing ( |10|) let us implement the Lagrange multiplier constraint 
so that Trn 2 = 1. Then 



TrnAn = 0. 



IF 



Let Tij k = ^/heijk where is the Levi-Civita symbol (satisfying ei 2 3 = +1). 
Then the volume form is = yhcPx = Tij^dx % dx^ dx k /3! and the Hodge star 
operation is given in terms of t^. For example = t^F^ jl for a two form 

F. Upon using the Bianchi identity, the a-model constraint and (|ll]) one may 
verify that 



VhTrh ij 



1 



r m e u F kL ± (nd iU + An) 



-T jmn e u F mn ± (nd jU + Dj 



n 



Vh Qe 2u Tr FijF ij + Tr D t nD l n + V*uV,-uJ ± 9, Tr {Vhe u T ijk F 3 



k n 



'12) 



A comparison of (|T2[) with (|S|) shows that we may rewrite the reduced action in 
the form fllPf) upon setting 

V(x\x 2 ,x 3 ) = e u , (13) 
and choosing units so that AnG = 1. We have defined vf in fllPP by 

= ^e u F fe/ ± (nd iW + An), (14) 

while the surface term B ± is 

B± = \J d ' ixd i (Vhh ij d jU ± Tr (VfreV^F^n)) . (15) 

Equations (0, 14) may be succinctly written in terms of differential forms as 



S 



]R {3 \h)r h - -Trv^ A *w ± + id(*dw ± Tr(2e u Fn)) 



(16) 
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with 



v ± = Vi dx l = *e u F± (ndu + Dn). (17) 

We remark that the requirement AirG = 1 sets a mass scale for the model. It 
has the same effect as in classical Newtonian mechanics of allowing static con- 
figurations of self-gravitating point charges when the masses and charges agree 
in suitable units. It is analogous to the choice of a critical coupling when con- 
structing Bogomolny equations for the Abelian-Higgs model of superconductors. 

As we commented upon earlier, the quadratic appearance of v ± in the action 
means that this term has an extremal for any field variation when = 
(for every % and gauge index a). These are the desired Bogomolny equations. 
In particular, supposing = 0, then varying the action with respect to the 
three-metric hij yields, 

5 J d 3 xVhR( 3 \h) = 0. 
This shows is a solution of the three dimensional Einstein equations 



R ( §\h)- 1 -h tJ R^(h) = 0, 



which in turn implies 



R^(h) = 0, R^\h) = 0. 



V) 



Using the three dimensional identity 

RijH — hikRji — huRjk + hjiRik — hjkRn + - (huhjk — hikhji)R 

we deduce that the three-metric is necessarily flat. This shows the various 
flat models assumed by previous authors were in fact necessary. If space is 
topologically R 3 then this flatness means it is isometric to Euclidean space and 
we may take = 5ij, in which case the four-metric is 

ds 2 = -V 2 dt 2 + A^x • dx. (18) 
V 

Let us summarise our calculations thus far. Assuming n satisfies Trn 2 = 1 
(and with a choice of units such that AnG = 1) we find 

n a d lU + D, t n a = T ± Tm e u F aM , (19) 
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or equivalently 



ndu + Dn = =p * e u F, 



(20) 



yield solutions of the field equations if and only if the three-metric hy 

defined by (|7|) is flat. These are our Bogomolny equations for the action (||); 
one can readily begin with ([II]) or (|20|) and a flat three-metric and derive (§-§). 
We have then established our first result and shown the compatibility of the 
Bogomolny equations with the remaining Einstein equations, so justifying the 
ansatz of the various authors mentioned in the introduction. 

As we mentioned in the introduction, the spacetimes with metric (|i~8|) are pre- 
cisely those static space-times for which super-covariantly constant spinors can 
exist ||13|| , and Tod has extended this classification to stationary metrics. 
Hlousek and Spector [|TJ]] have advanced arguments associating supersymmetric 
extensions to theories exhibiting Bogomolny bounds. This suggests a super- 
symmetric extension exists to the theory here under consideration. Here we will 
neither pursue the construction of such a theory nor address the issue of whether 



super-covariantly constant spinors exist; the recent work of |16[ looks at various 
4 dimensional supergravity reductions to sigma models addressing related issues. 



3 Solutions 



It remains to discuss solutions of the Bogomolny equations (19, ^). Before 
focusing on the case of axially symmetric solutions for the gauge group SU(2) 
we first show that an auxiliary magnetostatic problem Q may be associated with 
the equations and then discuss the case of covariantly constant solutions. 



3.1 General Properties 

First then let us consider projecting the Bogomolny equations in the direction of 
the cr-model field. This is analogous to projecting the usual Bogomolny equations 
of Yang-Mills-Higgs theory in the direction of the Higgs fields. Using (|i~T|) we 
find 

d l {e- u ) = ± l -r m Ti{F kl n), (21) 

2 Equally one may consider this to be an auxiliary electrostatic problem. Because we have 
in mind the nonabelian problem where the sources are magnetic monopoles, we adopt the 
magnetostatic perspective. 
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(equivalently de u — ± * Tr(Fn)) and this together with the Bianchi identities 
yields 

V k V k e~ u = ±~r ijk Tr (F i:j D k n). (22) 

Upon setting bi = di(e~ u ) and b = (bi,b 2 , b 3 ) we may recast our equations in the 
form 



V • b = p, V x 6 = 0, V 2 $ = p. (23) 



p = ±-T ijk Tr(F ij D k n), 6 = V$, $ = e~ u . (24) 



Here 



We see then from fl23|) that there is an Abelian magnetostatic problem associated 
with (|T^) in which $, b, and p play the role of the magnetic scalar potential, the 
magnetic field and the magnetic charge density, respectively. In particular the 
original metric coefficient goo is entirely determined by the Poisson equation of 

We observe that for the special sub-class of covariantly constant scalar solu- 
tions, those characterised by 

D(a a = 0, (25) 

the projected Bogomolny equations ( [2i~D in fact encode the full Bogomolny equa- 
tions. In this situation the magnetic charge density p (the source in the Poisson 
equation above) vanishes and the metric coefficient g 00 is determined by the solu- 
tion of a Laplace equation. Further, b is divergenceless in this particular situation 
and so it is the curl of an Abelian vector potential which has string singularities 
defined in the associated three-dimensional Euclidean space. The general solu- 
tion for V in this situation describes a system of point-like monopoles in static 
equilibrium interacting with the gravitational field. These are the Majumdar- 
Papapetrou solutions || 0. With only one monopole placed at the origin of 
coordinates this (positive mass) solution takes the form 

e " = 77m- (26 > 



where Q is the magnetic charge and r = ^/(x 1 ) 2 + (x 2 ) 2 + (x 3 ) 2 . The asymptotic 
form of the metric is given by 

Iftol « 1 " — , (27) 
r 

with a total gravitational mass M = Air\Q\. Such a solution corresponds to an 
extreme Reissner-Nordstrom black hole in which the null hypersurface r = 
corresponds to the horizon with radius \Q\ in Schwarzschild coordinates. 
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3.2 su{2) Solutions 

In order to proceed further with our analysis of the Bogomolny equations fllTf ) we 
restrict our attention to the su(2) case. When the scalar field is non-covariantly 
constant the full nonlinearity of the Bogomolny equations becomes apparent. In 
order to solve these equations we start by constructing a static, axially sym- 
metric, magnetic ansatz for the components of the Yang-Mills connection. Upon 
combining the resulting first order equations we shall obtain a single second order 
partial differential equation ( pE2"| ) governing the system and from whose solutions 
the various fields may be reconstructed. The restriction to spherical symmetry 
reproduces the ordinary differential equation obtained by earlier authors. Un- 



fortunately no analytic solutions for fl42|) (or indeed its spherically symmetry 



reduction) are known and we must consider various approximations. 

In the su(2) setting a non-covariantly constant Din a may be expressed as 

D in a = a lP a + (3 iq \ (28) 

where ati and $ are functions of the space coordinates and p a and q a are cho- 



sen to satisfy p a n a = q a n a = 0, so that ( |TT| ) holds. Imposing the additional 
conditions p a p a = q a q a = 1 and q a = e abc n b p c , the triad (n,p,q) becomes a ro- 
tating, orthonormal base for su(2). Working in spherical coordinates, we adopt 
an ansatz in which we assume that n, p, and q take the special forms 

n = sin 9 cos T\ + sin 9 sin T 2 + cos 9 T3 = x a T a , 

p = cos 9 cos T\ + cos 9 sin T 2 — sin 9 T 3 , (29) 

q = — sin Ti + cos T 2 , 

where Ti,T 2 ,Ts are the generators of the adjoint representation of s«(2).[] 

Now upon combining the definition of the covariant derivative Din a with 
(2"%p we may solve (extending pUfl) for the Yang-Mills potential A — A4 dx l = 



A? T a dx l , to obtain 

A a = a t q a - (3 iP a - e\ c n b d in c - <5,n a , (30) 

where Si are three arbitrary functions of x l . This means the components of the 
gauge field may be expressed in polar coordinates as 

A a = q a {a- d9) +p a {sin 9d<f) -(3)- n a 5. (31) 



3 With [T ai T b ] = e abc T c we have 

q=[n,p], n=[p,q], p=[q,n], 

and 

dn = pd9 + sin 9q dtp, dp — ~nd6 + cos 9q d(f>, dq = — (sin On + cos 6p)d<p. 
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By projecting ( |20| ) in the (n,p,q) directions we find the equations connecting 
these various unknowns (a = cii dx l , etc.) to be 



du = =F * e u (—d5 + a A (3 — sin 6d6 A defy , 
a = =F* e u (—dp + 5 A a + a A cos6d<p) , (32) 
/3 = =F * e M (da — P A 5 + P A cos ( 



Further, we have that 



[L 3 - T 3 , A] = q d^a - p d^P - n 8^8, (33) 

where L 3 is the generator of space- rotations around the x 3 axis. By taking a, P 
and 5 independent of <fi we then have 

[L 3 -T 3 ,A]=0. (34) 



This property is characteristic of an axially symmetric, su(2) connection ||2T . 



At this stage we have expressed the Bogomolny equations (|i9D for the gauge 
group su(2) in the form fl3"2]). No approximations have so far been introduced 
and the problem remains of either solving (|32|) or finding ansatz that enable their 
solution. 

Our ansatz for the study of ([19]) is now based on the following choice for a«, 
Pi, and Sf. 

ai = 0, a 2 = x, «3 = 0, 

Pt = 0, p 2 = 0, p 3 = sin6 X , (35) 
6 1 = , 5 2 = , 5 3 = ij, 



where x an d ^ are m general functions of r and 9. These choices mean ( |34| ) 
is satisfied and we have an axially symmetric, su(2) connection. Spherically 
symmetric solutions are covered within our ansatz by ip = 0, x = x( r )i m which 
case we recover the t'Hooft ansatz, 

A1 = e% h -(l-x(r))- 

r 



4 It is perhaps worth remarking that an axially symmetric connection does not necessitate 
an axially symmetric space-time. This will occur when V (or equivalently e u ) is independent 



of 4>. From (32) we see this means 

= d r 5 2 - doSi - ai/3 2 + a 2 Pi, 

and so for (Q) to admit d/dcj> as a Killing vector implies a constraint on some of the connection 
parameters. Similarly further constraints are needed for the metric to be spherically symmetric. 
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With this ansatz and taking^ the (-) sign in ( |T9| ) we find five distinct equations 



amongst the components of the Bogomolny equations (^2|): 

d 



()r X = -e~ u X (36) 



d ipx 

^ = ^"H: ( 37 ) 

oO sin 6* 



g r ( ) r 2 ( ^ sin^ (96 1 / ' 



' sintfcV 1 J 



I: ( e "") = °- ( 4 °) 



Not all of these equations are independent: ( p9|) is a consequence of ( p6|) and 
(|37|). Now equation (pEPP implies that our conformstatic metric ([/]) admits ^ as 



a Killing vector. Within our ansatz the a-model equations (§) reduce to (|3 
and the axisymmetric charge density p ([24|) simplifies to 

1 d(v 2 ) 

p = 7*fr- (41) 

The latter result indicates that e~ u = y is not in general an harmonic function. 
Combining yields the equation 

r2 ^ lnW + ^( 8i,1 ^" lW ) = 1_x2 ' (42) 

This equation governs the solutions of our Bogomolny equations. Once we have 
a solution for (fj2|), the remaining unknowns V and ip can be simply determined 
using (pD, (HD and dp. 

At this stage we have obtained the single second order partial differential 
equation fl42|) governing the system. In the case of spherical symmetry this 
equation reduces to a form of Emden's isothermal gas equation. (In this case 

5 Choosing the opposite sign leads to the same class of solutions. The sign ambiguity can be 
absorbed by working with a quantity ±e M ; this sign is responsible for the existence of solutions 
with either positive or negative magnetic charge and positive gravitational mass. 
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the a-model equations (|5|) become trivial as (|3"7|) vanishes.) With ln|x(r)| = 
—F(r) — In r we obtain equation (21) of 

r A F" = e~ 2F . (43) 

Such Emden type equations may be rewritten in the form of an Abel equation 

dx rx 



dr r — 1 + x 



2 ■ 



(44) 



where r = r F'(r) — 1 = r e " (and x = e F l r )- Indeed the variable r naturally 
arising in this transformation is such that the metric takes the form 

' 7 \ 2 

,2 2u t.2 , -2u I ar \ 7-2 , -2/7/)2 , -2 a J a2\ 



ds 2 = -e 2u dt 2 + e~ 2u 1^1 df 2 + f 2 (d6 2 + sin 2 6d<f) 2 ). (45) 

Unfortunately no analytic solutions for (fj2]) (or indeed the spherically symmetry 
reductions ( |43"D or ([HQ) are known and we must consider various approximations. 
The form of the metric (|45D is helpful when considering possible singularities of 
solutions, with 47rf 2 being the area of a two-sphere about the centre of symmetry. 
We see from this that the conditions for a horizon are 



du 
dr 



= 1, 

H 



= 0, (46) 

H 



and consequently r = at a horizon. 

It follows from our equations that an asymptotically flat, spherically sym- 
metric black hole solution with regular gauge fields is necessarily an extreme 
Reissner- Nordstrom black hole. This may be argued as follows.^ The spherically 



symmetric Bogomolny equations together with (46) entail x\ H = f° r regular 
gauge fields. Further, the asymptotic boundary condition on the gauge field is 
x\ 00 = 0- Now it follows from (^) that x is a monotone function from the 
horizon to radial infinity, and so x{ r ) vanishes identically. This brings us to 
the covariantly constant extreme Reissner-Nordstrom solution of the previous 



section. This result agrees with the work of |23[ who have analysed the region of 
parameters for which spherically symmetric EYMH black holes exist: our sigma- 
model corresponds to their (/3 = oo) infinite mass limit, while the Bogomolny 
equations necessitate their remaining parameter to be a = 1. Without the re- 
striction of spherical symmetry it is more difficult to determine the singularity 
structure of our solutions. The Majumdar-Papapetrou solutions of the previous 
section may be shown to have singularities shielded by horizons. We will defer 



3 This argument is based on one shown to us by E.J. Weinberg and D. Maison. 



12 



such discussion on the nature of singularities until the appendix, looking here to 
the possible asymptotics of solutions. 

If is small for large r, then fl42| ) and (^) have the approximate (asymp- 
totically flat) solution 

x (r)«B2?, e««l-± (47) 



where B is an arbitrary constant. Comparing this asymptotic solution with (|27| ) 
we conclude that ( |47D corresponds to a magnetic monopole with unit magnetic 
charge and total gravitational mass M = Air. Another approximate solution for 

(ED ^ 

In \x\ ~ Cv^ 7 sin I — lnr + Q j , (48) 

where C and Q are integration constants. As discussed by Balakrishna and Wali 
||, fl4"8] ) is valid for \x\ ~ 1, which is the case when r — > 0. The oscillatory 
behaviour of ( |48| ) implies a countable, infinite set of (naked) singular spheres 
surrounding the origin, which is the locus of an event horizon of zero surface 
area. (See Q for details.) We have found numerical evidence for the existence 
of asymptotically flat solutions of (^) which smoothly interpolate between the 
(asymptotic) monotone and oscillatory regimes described approximately by (|47]) 
and (f48"D, a point left unanswered in earlier works. This is described in the 
appendix alongside remarks on the nature of the singularities associated to such 
solutions. 

In the non-spherically symmetric case, we may begin with the approximation 
(for large r) to (f42|), 



e ' 

X « B — 

r 



1=1 



(49) 



Here B is is arbitrary and e is a small parameter. Upon using (|36|) and dropping 
terms of order e 2 we find 



1 _ + e ylA Pll 

iy> ^^^^ cr>\-\~\ 



^COSl 



(50) 



Comparing with fllTf ) we see this corresponds to an asymptotically flat solution. 
Similarly using (|37D we obtain 

. Q ^ fi dPi(cos6) 
tfia-eBmO^-j ^ . (51) 

i=i 
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Observe that, in contrast to x> ip is a long range potential. It is also possible to 
find an approximate solution of ( 42|) for small r, such that \ is finite at r = 0: 




\\ wl + CVrsin — lnr + fi + V q r 5+VW +1 K p,( C os 0), (52) 



where C, f2, and Ci, C2, C3, .... are arbitrary constants. This expression shows a 
remarkable difference between the radial dependence of the monopole and higher 
multipole terms in \- We see that fl4T3| - |5"2D reduce to the previous spherically- 
symmetric results when every /; and q vanishes. 

A simple perturbative method can be used in order to improve the asymp- 
totic solutions ([47|) and ( f4l^|51| ). We illustrate the procedure in the spherically 



symmetric case only, but it can be extended straightforwardly to the case in 
which x depends on 9 as well. We have for large r that x{ r ) ~ Be~ r /r which 
gives the correct asymptotic behaviour for the metric. Now suppose that 

— T 

X (r) = B—(l+e) 
r 



where £ C 1. Upon substituting in fl42|) we find the approximate expression 

B 2 e- 2r 



(53) 



X*B— (1-- — ), (54) 



4 r 4 

As a consequence, our new asymptotic solutions are 

.e^ / & e~ 2r 
r \ 4 T 



2m(r) , , 

^00 « 1 55 

r 

where 

B 2 e~ 2r 

m ( r ) = 1 - -r, — r- ( 56 ) 

2 r 6 

This improved asymptotic approximation exhibits a very small correction in 
the distribution of gravitational mass, suggesting the existence of a massive, 
extended magnetic core. Such exponential corrections also arise from instanton 
corrections. 
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4 Inclusion of a Dilaton 



We now extend our earlier analysis to include a dilaton. Up to a redefinition of u 
we will find that the Bogomolny equations are unaltered. Consider the EYMDcr 
system given by the action 



S, 



d A x\/\g\ 



16nG 



- Tr {D^nD^n) + - (Tr n 2 - l) 



(57) 



The field equations are now (0) and @ together with the modified gauge equa- 
tion 



(5f 



V\9\ V J 
and the new dilaton field equation 

j=d, (y\9\d^) = (47rG7) e lv Tr F fMU F tMU . 



9\ 

Upon making use of (|S|) one observes [|nj that with the identification 

<P = 2 7 In V\9oo\ 



(59) 



(60) 



the dilaton equation of motion ( |59|) coincides with the (OO)-Einstein equation. 
Our construction of a reduced action now proceeds as before. To employ the 
identity (0) we want 



„2m 



i.e. e" 



which, upon using (|60|), requires 

u = (l + 7 2 ) InVbool 
Finally, the choice 

47rG(l + 7 2 ) = 1 
means we have a reduced action 

d 3 xVh 



^LRW(h) --Trti j v?vf 
4 y 1 2 13 



+ B, 



j ■ 



(61) 



(62) 



(63) 



where is again a surface term and vf is given by (|I~4"D. Up to the scaling of 
the scalar curvature this is of an identical form to (|TUD. In particular this means 
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we obtain the same Bogomolny equations as previously with the same conclusion 
that the three- metric is flat. 

Therefore, assuming n satisfies Trn 2 = 1 and with a choice of units (§2), we 



find that the Bogomolny equations ( |19| - |2~0| ) again provide solutions to the field 
equations for fl57D if and only if the three- metric hij is flat. The dilaton field is 
given by (|59|) for such solutions. 

An important consequence of having the same Bogomolny equations is that 
the solutions discussed in the previous section may be directly used in the present 
setting, though with the metric function V now related to e u via (^) rather than 
(P~5|) . Thus in the absence of a dilaton the solution to Laplace's equation 
describing a single monopole led to the asymptotic form of the metric (p7| ) we 
now find 



000 



1 - 



2m. 



m 



\Q\ 



Whereas without the dilaton we have an extreme Reissner- Nordstrom black hole, 
according to Ref. |J the solution with dilaton has a naked, point-like singularity 
at r = for 7 ^ 0. The energy integral may be determined analytically for this 
solution with arbitrary 7 yielding 



E 



27r(l + 2 7 2 )|Q| 
I+7 2 



The fact that E is finite in the case 7 7^ -in which the integration region includes 
the naked singularity- can be considered as the extension to curved spacetime of 
a result previously found by Bizon |l7j , who obtained a Bogomolny- type solution 
for the Yang-Mills-dilaton equations. (See also [IS] for an independent discussion 
of these equations.) Park [|19|] has also discussed the naked singularity structure 
of dilatonic point-like Taub-NUT multi-body solutions to an Einstein-Maxwell- 
Dilaton theory. 

Similarly with our axially symmetric solutions the metric function V is mod- 
ified accordingly, based on the same function e u . Thus to a first approximation 
|) leads to 

1 1 



V 



which may be improved via (55) with 



m(r) 



1 + 7 2 r 



B 2 e 



-2r" 



1 + 7 2 



16 



5 Solutions with Higher Winding Numbers 



In this Section we study axisymmetric solutions of the modified Bogomolny 
equations with winding number greater than one. Again the problem lies in 
finding tractable ansatz for the exact Bogomolny equations. To this end, let 
us consider the following [I 



21fl magnetic, static prescription for the SU(2) 



potentials: 



dx" = q {k) 



-—dr + (H 2 -l)d9 



+ k [n {k) H 3 +p {k) (I-H4)] sin 



where n k , p k and q k are defined by 



n ( k ) = s j n q cos ^ 7^ _|_ s j n Q s i n t 2 _|_ cos Q T 3 , 
pW — cos cos kcp Ti + cos 9 sin k<p T 2 — sin 9 T 3 , 
= — sin k<p T\ + cos kcj) T 2 , 



(64) 



(65) 



and the four functions Hi depend on r and 9 only. The integer k will be inter- 
preted as the winding number of the solutions. (Equations ([64]) and ( |65"D have 
been inspired by the ansatz considered in [21], which is set up in the fundamen- 
tal representation of SU(2).) This new prescription for the connection can also 
be obtained as a consequence of (|30D if we replace^] n, p, q by rv- k \ p™, q^ k \ 
respectively, and make the following choice of parameters: 



Qti = — 

Si = 0, 



Hi 



a 2 

fa 

5 2 -- 



■H 2 , 
= 0, 
0. 



a 3 = 0, 

03 = k sin 9 if 4 
5 3 = —ksm9 H 3 . 



(66) 



This ansatz satisfies the axial symmetry condition 

[L3-kT 3 ,Ai] = 0, 



which generalises ( |34] ) to higher winding numbers. We also observe that, for 
k — 1 and Hi — 0, H 2 — H4 — x(r, 9), H 3 = — ^pj- , our original axially 
symmetric ansatz (p5|) is recovered. 

Combining this new ansatz for the connection with the projections of the 
Bogomolny equations on the rotating basis (|65|), we obtain the following system 
of equations: 



r 



r)H 

- +cot# (H A -H 2 )-H 2 H 3 



09 



(67) 



Now dn^ 



o^de + k smOqWdtj), dpW 



-k (sin On^ + cos 9p^)cbf> and A = qW (a - dO) + (k sin 



-n w d9 + k cos 6><? (fe) dg (fc) 



/?) - n ( fc )<J. 
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H 2 e 



-u Tk 



r—^ + cot8H l + H X H 3 
or 



(6£ 



dr r d6 



(69) 



9r 



n tt 

o9 



(70) 



ar r 



(71) 



90 



(e~ u ) = 0. 



(72) 



An asymptotically flat, covariantly constant solution of ( 167[ - |7"2"|) can be easily 



obtained if we assume H 1 = H 2 = = and H 3 
function e~ u is given by 



sin 9 



. The harmonic 



\Q[ 

r 



'MCOSI 



where Q = k is the total magnetic charge. The corresponding solution for ip is 

(73) 



, , sin# ^-^ Ci dPiicost 

^ = ^-—1^1?- 

i=i 



d0 



The identification of total magnetic charge with the winding number in this 
solution implies that the gravitational mass is conserved for topological reasons. 
(See [0 for a discussion of this point in the case with no dilaton present.) We 
have yet to solve this ansatz for non-covariantly constant cr-model solutions. 



6 Discussion 

This paper has examined the Bogomolny equations and their solutions for Einstein- 
Yang-Mills-a models (with possible dilaton couplings) on static space-times. Our 
derivation of the Bogomolny equations leads to several new observations. In par- 
ticular the Bogomolny equations are consistent with the Einstein equations if and 
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only if the associated (conformally scaled) three- metric is flat. These are pre- 
cisely the static metrics for which super-covariantly constant spinors exists, and 
this class of metrics includes all of the particular ansatz considered by previous 
authors. (Stationary metrics for which super-covariantly constant spinors exist 
have also been classified, and the extension of such models to these space-times 
will be dealt with elsewhere.) The connection with super-covariantly constant 
spinors suggests a supersymmetric extension of these models that we have not 
so far pursued. Although we have assumed a cr-model interaction throughout 
we do note that recent results [^5| for 2+1 gravity coupled with non-Abelian 



Chern-Simons vortices on an axially-symmetric spacetime indicate Bogomolny 
equations are possible (subject to appropriate ansatz) for particular (and in this 
case nonrenormalizable) potentials. It is an interesting question to see whether 
the methods of this paper will allow an extension to nontrivial potentials. 

Having obtained the Bogomolny equations we have also considered their so- 
lutions. Two cases arise depending on whether the cr-model field is covariantly 
constant or not. The former situation is somewhat easier to analyse and leads 
to a set of (Euclidean) Abelian magnetostatic equations, valid for any gauge 
group. These equations are sufficient to determine solutions of the Bogomolny 
equations in this case. 

When the cr-model field is not covariantly constant the Bogomolny equations 
are rather more complicated. Here we have focussed on axially symmetric solu- 
tions for the case of su(2), though the situation involving larger gauge groups 
is very interesting. By projecting onto the various gauge components we may 
express the Bogomolny as the coupled system of first-order equations fl32|) . The 
imposition of axial symmetry on space places additional constraints on the pa- 
rameters of an axially symmetric gauge connection. At this stage we adopted an 
ansatz^] to simplify the analysis of (p^) and we recovered the results of previous 
authors who studied the spherically symmetric case. In fact we were able to show 
that an asymptotically flat, spherically symmetric black hole solution with reg- 
ular gauge fields is necessarily an extreme Reissner-Nordstrom black hole. Our 
work also shows a totally different type of radial dependence between the higher 
multipoles and the monopole terms in our axially symmetric solution. Inter- 
estingly an improved asymptotic approximation suggested a 'massive magnetic 
core', the asymptotic corrections having the same form as instanton corrections. 
In our final section we extended the original ansatz to include higher winding 
numbers finding an exact, covariantly constant solution for arbitrary winding 
number. We have not succeeded thus far in finding solutions to this ansatz 
for non-covariantly constant solutions with higher winding numbers. Although 



8 One hope is to find an integrable system in these equations paralleling the appearance of 
the Toda equations in the flat space monopole equations. 
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obtaining exact, nonlinear solutions may be an impossible task, expansions in 
powers of - and 6 (in the spirit of f24|) could help to understand some aspects 
of the solutions. Such an approach is left for future work. 

Clearly obtaining the Bogomolny equations is far from the end of the matter. 
Although this system of equations is of lower order and simpler than the original 
field equations, they are still complicated. Our analysis has presented some of 
the features of these equations and we look forward to their wider study. 
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Appendix: Remarks on Singularities 

In this appendix we consider the possibility of finding non-singular solutions to 
the gravitationally coupled Bogomolny equations, restricting our attention to 
the case of vanishing dilaton (7 = 0). 

Besides the invariants R, R a pR al3 ', and R a ^^ v R v , the simplest one that can 
be used in the analysis of singularities is J = F^ v F a ^ v , which is the non-Abelian 
generalisation of the one used by Hartle and Hawking in their study of Abelian, 
multi-black hole solutions Using ( |i9|) we obtain the expression 

J = 2 (D i n a D i n a + dtu&u) . (74) 



If we assume that n is covariantly constant and choose the metric fllq) , then (|74| ) 
reduces to J = 2 (VU) 2 . In this special case 1/V is harmonic (see Section 3) 
and can be determined exactly even without spatial symmetries. One solution 
is 

p=i k 

where p p = \J [x — x p ) 2 + (y — y p ) 2 + (z — z p ) 2 . This solution describes a collec- 
tion of N point-like monopoles with magnetic charge m p , located at (x p ,y p , z p ). 
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For stationarity every m v has the same sign and ([75]) describes the Majumdar- 



Papapetrou solutions considered in |2(J and J is finite everywhere except at 
some points which are enclosed by event horizons. Also, according to Hartle and 
Hawking, other distributions of charge different to discrete point sources can be 
considered. However, the corresponding solutions can be shown to have naked 
singularities. We conclude that the covariantly constant solutions of ([19]) are 
necessarily singular, but these (genuine) singularities are not naked when the 
source is a discrete set of (like charged) point-like magnetic charges as described 
byfl). 

The analysis of singularities in the non-covariantly constant case is much 
more difficult because 1/V is not harmonic and we do not know any exact, non- 
trivial solution of the non- linear equations (|T9"|). Nevertheless, we present some 
considerations based on our knowledge of asymptotic solutions. 

Restricting our attention to the axially symmetric case, we use spherical 
coordinates and (|35|) to obtain the reduced expression 



J = JLJL + 2 (W) 2 . (76) 

Both terms in (|76|) are non-negative, so there is no possible cancellation of in- 
finities and, exactly as in the covariantly constant case, the solution is singular 
wherever 1/V = 0. 

Using (PS|) and flTB|) it is easy to see that J is unbounded at each 



one of the infinite critical points of ln|x(r)| near r = thus yielding genuine 
singularities. One may ask whether there exist spherically symmetric solutions 
of (|42|) with such pathological behaviour for very small r yet that behave well in 
the far asymptotic region. We have numerical evidence for the existence of such 
a class of spherically symmetric, asymptotically flat solutions. The transition 
between the (asymptotic) monotone and oscillatory regimes of ( [47|) and (|48| ) 
can be conveniently analysed in terms of the new variable £(x) = ln% 2 , where 
x = lnr, which satisfies the equation 

^-§- + 2e* = 2. (77) 

ox z ox 

This transformation enables a simple numerical analysis near x = 0. The solution 
shown in Figure 1, satisfies £(0) = 1, £'(0) = 0. It evidences the different 
behaviours of fl47|) and (|48]), and displays a smooth regime transition at x — 
(r = 1). As a consequence of fl36|), V 2 is infinite at every x where £'(x) = 0. It 
appears then that (|42|) admits at least one class of asymptotically flat, non-linear 
solutions which contain essential singularities. 

We next would like to extend our analysis of singularities to non-linear so- 
lutions in the non-spherically symmetric case. This task is even more difficult 
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because at present we only know the asymptotic behaviour of putative solutions 
and have not answered the question of existence. Equations fl49D-(|5lD show the 
compatibility of fl4"2] ) with an asymptotically flat spacetime, and (^) exhibits 
the asymptotic behaviour of any possible solution satisfying |x( r ;$)| ~ 1 near 
r = 0. As in the spherically symmetric case, the (oscillatory) monopole term in 
ln|x(r, 9)\ leads to spacetime singularities which can be avoided only if C — 
in (|52|) . Also, it is easy to see that the new leading I = 1 term makes V behave 
as y^gj when r — > 0. This asymptotic behaviour also implies unbounded values 
for J near r = as a consequence of (|76|) , so we also put c\ = in ( ^2|) and 
proceed with the analysis of the remaining terms in that series. The asymptotic 
behaviour of the higher multipole terms (7 = 2,3,4, ...) once more implies that 
1/V vanishes at r = 0, and so at that point J is infinite again. Clearly, this 
implies the existence of spacetime singularities arising from solutions which have 
the asymptotic behaviour |x(r, 9)\ « 1 near r = 0. However, we can say noth- 
ing about the existence of such singularities in solutions with a different type of 
asymptotic behaviour. We suspect that every asymptotically flat solution must 
give rise to spacetime singularities near r = 0. We doubt that analytic methods 
will be useful at this stage for studying the solution space of (j^), so the next 
step should be a detailed numerical study of possible non-spherically symmetric 
solutions. 
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Figure 1: Numerical solution for equation 
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